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Abstract 

Discussed is the structure of classical and quantum excitations of in- 
ternal degrees of freedom of multiparticle objects like molecules, fuUerens, 
atomic nuclei, etc. Basing on some invariance properties under the action 
of isometric and afBne transformations we reviewed some new models of 
the mutual interaction between rotational and deformative degrees of free- 
dom. Our methodology and some results may be useful in the theory of 
Raman scattering and nuclear radiation. 

Introduction 

Roughly speaking, in what follows we are dealing with some ideas concerning 
applications of affine symmetry in fundamental physics. The very idea is not 
new. Let us mention, e.g., some papers by Hehl, Ne'eman, Sijacki and others 
(see [Ill3nil3Il[SIllIl[7i|771[7SllSl]and references therein). 

Everybody knows from ancient times that geometry of the physical world 
is organised in a hierarchical way and consists of some modules. The most el- 
ementary of them is one based on the ideas of Thales of Miletus. The crucial 
concepts and ideas there were just the parallelism and proportion of segments 
on the same straight line or, to be more precise, on parallel straight lines. Much 
of our geometry may be developed in this way. And further on, the "brutal" 
metrical concepts of distances and angles, or equivalently — scalar products, 
did appear. Mathematically speaking, they have to do with bilinear forms on 
the space of translations in the physical space. No doubt, the metrical concepts 
are very important also on the fundamental level. Indeed the very pseudo- 
Euclidean metric structure of the Minkowskian four-dimensional space-time or 
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the Euclidean three-dimensional space are, so to speak, a vacuum-state mani- 
festation ("ground state" in a sense) of the very important field, i.e., just the 
gravitational field, the force known to humanity from ancient times. Neverthe- 
less, if the Thales-Euclid hierarchy of concepts is seriously treated, it seems very 
natural to ask what would be the corresponding hierarchy of physical concepts 
and theories. Then one can ask, even if motivated by a purely academic in- 
terest, what would be a hypothetical afline physics and what natural schemes 
of the metrical (isometry-based) breaking of that Thales-preestablished affine 
harmony seem to be physically and experimentally justified. And why at all? 
Is this some additional interaction or something like the spontaneous symmetry 
breaking? 

The idea of using GL(3,M), SL(3,M) or other non-compact groups, e.g., 
symplectic and pseudo-orthogonal ones, is not new in fundamental physics, let 
us refer again to Hehl, Ne'eman and Sijacki et. al. [131 [30l [3ll [HI [521 [76l [711 [78] . 
The pseudo-orthogonal group O (2, 4) is related in a known way to the theory of 
hydrogen atom. In nuclear physics the group GL (3, K) (or rather its subgroup 
SL(3,M); nuclear matter seems to be as incompressible as usual macroscopic 
fluids) has been used for a long time in the collective droplet model of nuclei 
|ni[ini[ISl[I3[5Tl[ini[HDl[lIl[Sl[Ml[Ml[HIl[HSl- Some of the mentioned non- 
compact groups were thought on as configuration spaces. They were also used 
as so-called non-invariance groups, spectrum-generating groups (or algebras), 
e.g., Hamiltonian was one of the generators, and the spectrum might have been 
found in a rather algebraic way, on the basis of commutation relations alone, 
etc. (see, e.g., [3 [3 [39l [50] and references therein). 

In a sense the classical model of affine degrees of freedom has been used 
for ages in the theory of the Earth's shape (Riemann, Dedekind, Maclaurin; it 
seems that in a sense even Newton himself; see, e.g., [Tnil3Hl[Zl] and references 
therein). In XX-th century the interest in this topic revived, as mentioned, in the 
connection with nuclear and molecular dynamics [42l[66l[67l[72l[73], astrophysics 
[T5] and the theory of continua with microstructure [Ti l H5 1 [ii l [72 1 [75 1 [57] . 
The famous concept of ellipsoidal figures of equilibrium became fashionable. 
On the purely mathematical level this effort was accompanied by new methods 
based on Poisson geometry and the dynamics on co-adjoint orbits of Lie groups. 
The subject has been intensively studied both on the classical and quantum 
levels [a [3 [Ml [351 [Ml [Sana [571 [721 [7g. 

As mentioned, we concentrate on the linear group as one which rules de- 
grees of freedom of internal/relative motion. So, the only admissible modes of 
motion are spatial translations, rotations and homogeneous deformations. In a 
sense, our topic belongs to the dynamics on Lie groups as developed by Arnold, 
Marsden, Hermann and others [H [31 [311 [23] , in particular, something like "afline 
Euler equations" appears [721 [73] . But it is so " in a sense" only. Degrees of 
freedom are ruled by Lie groups GL(3,R), GAf (3,R), but it is not so with 
Lagrangians, metrics on the configuration space and equations of motion. But 
the very taste of systems with group-theoretic degrees of freedom consists in 
considering left or right (or both, i.e., two-sides) invariant geodetic models. It 
happens quite often that some explicit solutions may be analytically obtained 
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and expressed in terms of exponential mappings on Lie algebras or in terms of 
some well- investigated special functions on groups [34l |35l [72l [73l [83l [84l [85] . 

But for all applications and models mentioned above the afEne or linear 
groups do not preserve Lagrangians, Hamiltonians and equations of motion. 
The corresponding metrics on the configuration space have much weaker, at 
most isometric symmetry (Euclidean motions, preserving the Euclidean metrics 
of the physical space). So, and this is really disappointing, there is no profit 
of using group-theoretic degrees of freedom. Hamiltonian generators of the to- 
tal affine group in general fail to be constants of motion or good quantities to 
be used in balance laws when external influences are taken into account. And 
geodetic models, i.e., ones without forces and potentials, are not useful, because 
in non-compact spaces they, as a rule, predict the unbounded, "escaping to 
infinity", i.e., "dissociative", behaviour and no bounded motion, i.e., no oscilla- 
tory "elastic" phenomena. Similarly, on the quantum level one obtains purely 
continuous spectrum, no bound states. 

Let us compare all that with the beautiful structure of left-invariant geodetic 
models on the rotation group, i.e., with the translations- free rigid body (pure 
rotations) and with the right-invariant geodetic models on the group of volume- 
preserving diffeomorphisms, i.e., with the Arnold description of incompressible 
ideal fluids [1 H [3 [ni HOI HH iS] . And we try to repeat this reasoning for affine 
bodies. Quite independently of any perspective of physical applications, it would 
be interesting at least from the purely academic point of view of analytical and 
quantum mechanics to investigate the theory of systems with the kinetic energy 
(Riemann structure of the configuration space) just exactly invariant under spa- 
tial (Euler coordinates) or material (Lagrange coordinates) affine group. Let us 
also mention that the model right-invariant under SL (n, R) (under unimodu- 
lar and orientation-preserving transformations of Lagrange variables) may be 
considered as a very drastic discretization of the Arnold description of ideal in- 
compressible fluids. Discretization is drastic, because it reduces the continuous 
system with infinite number of degrees of freedom to one with a finite number 
of them, namely n(n + 1) in the n-dimensional "physical" space, i.e., 12 in the 
physical case when n — 3; among them 3 translational and 9 internal (relative 
motion) degrees of freedom (respectively n and in the academic n-dimensional 
world). Applications, e.g., to "nuclear fluids" seem quite possible. Affine invari- 
ance in Euler variables is a bit more doubtful and academic, nevertheless, there 
is some logical necessity which tells us that such models must be formulated 
(that is done below) and their consequences must be carefully derived and com- 
pared with the experimental data (the last thing yet to be done). In any case, 
quite independently on the academic idea of affinely-invariant models based on 
the Thales part of geometric axioms, there are some physical arguments which 
encourage us to undertake this effort. There are known situations where in 
condensed matter, especially in strongly condensed matter, the effective kinetic 
energy is not based on the "true" metric tensor, but on the tensor depending 
somehow on physical parameters. The most important example is the effective 
mass tensor of electrons in solid state physics [36] . The point is that the inertia 
of electrons becomes then less dependent on the spatial metric tensor and more 
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sensitive to surrounding physical factors influencing their motion [35]. Some 
phenomena like those are expected also in dynamics of defects in solids or in 
the theory of a bit exotic objects like bubbles and voids in continua [Tl], e.g., 
in fluids. The more so, such exotic fluids like nuclear matter may perhaps be- 
have in a way compatible with predictions of affinely left-invariant models (ones 
afhnely invariant in the physical space). In any case some interesting things 
are expected when, e.g., the effective mass tensor is proportional to the Cauchy 
deformation tensor. Geometrically this is the simplest model affinely invariant 
in the physical space. 

It must be strongly stressed, and the above introductory remarks prove 
this, that the models analysed below, in spite of their using of GL(3,R) or 
GAf (3,]R) and their subgroups, have httle, if anything, to do with otherwise 
very respectable and physically useful historical models discussed classically 
by Riemann, Dedekind, Dirichlet, Chandrasekhar, Bogoyavlenski (|15 ( \72 \ [73] 
and references therein), and on the quantum microphysical level by Weaver, 
Cusson, Biedenharn, Rosenstell, Rove and many others UHl HH IIZI HH HSl 
[501 HU EH ESI IHS]- Their dynamics are completely different. Because 

of this there are also some differences in analytical methods. When dealing 
with matrix groups, calculations are usually based on various decompositions 
of groups. For our purposes the most convenient ones are polar and, first of 
all, singular value decompositions; for latter one we often use the suggestive 
term " two-polar decomposition" . Other well-known multiplicative splittings 
like the Gauss or Iwasawa decompositions do not seem to be effective within 
the framework of our models. 

The peculiarity of our model is that, according to the Noether theorem, 
affine momenta, i.e., generators of affine group and their functions, are quantum 
constants of motion. More precisely, it is so for generators of the volume- 
preserving subgroup. In particular, if translational motion is neglected, the 
Casimirs of SL(n, R) are constants of motion [72l [73] . 

Another very interesting and important feature of affinely-invariant kinetic 
energies is that they may be used as Lagrangians describing the dynamics of 
large isochoric elastic vibrations coupled with rotational motion. No poten- 
tial term is necessary in such models and all advantages of invariant geodetic 
dynamical systems on Lie groups or their homogeneous spaces may be used 
[il3ll34ll35llMl|55l|56l|57l|72l|7g. The idea of encoding the dynamics in 
an appropriately chosen metric tensor of the configuration space resembles the 
Maupertuis variational principle, however, in our model the mentioned metrics 
are very natural and follow from some first principles based on geometric ideas. 
If the object is compressible, then dilatations must be stabilized by some po- 
tential term depending on the dilatational ratio, i.e., on detf in the formula ([1]) 
below. It is natural to assume something like the potential well or other model 
rapidly growing when detip deviates from unity [TH] . 

It must be stressed that our constituent point particles are spinless. The 
term "spin" in the classical and quantum parts of our treatment refers to the 
orbital angular momentum with respect to the centre of mass. When trying to 
apply the model in nuclear physics one must modify it by introducing in addition 
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and carefully taking into account the "true" spin of constituents, i.e., nucleons. 
Without such corrections our model may describe only some aspects of the 
nuclear dynamics. As a matter of fact, the same concerns molecular dynamics, 
although spin effects are then less essential than in nuclear phenomena. 

Incidentally, there are some subtle problems concerning the concept of spin, 
its afhne aspects, and the relationship between spin and orbital angular mo- 
mentum. And some ideas, perhaps rather speculative and hypothetical, may 
be formulated. More precisely, we mean here some speculations about the pos- 
sibility of "explaining" spin as the "orbital" angular momentum of internal 
motion of extended or composed objects. In standard theories spin is a primary 
characteristic of elementary particles, described formally in terms of irreducible 
unitary representations of the group SU(2) (Spin(n) in n dimensions) acting on 
multicomponent wave functions or fields subject to the quantization procedure. 
The use of the covering group of S0(3,R) (SO(n,]R) in n dimensions) results in 
appearing of half-integer spin, predicted on the purely experimental, spectro- 
scopic basis by Goudsmit and Uhlenbeck [5D] . This is angular momentum as an 
irreducible concept; angular momentum without rotation of anything. There 
were attempts of formulating afhne fundamental physics based on irreducible 
unitary representations of GL(3, M), the covering group of GL(3, M) (GL(n, R) in 
n dimensions) |521l76j. Of course, those representations are infinite-dimensional; 
the problem is additionally complicated by the fact that GL(n, M) is not a linear 
group (does not possess faithful realizations in terms of finite matrices) . When 
restricted to SU(2) (Spin(n)), those representations become reducible and split 
into direct sums of infinity of irreducible representations found by Wigner and 
Bargmann [30l[3Tl l52 l l76j: obviously, they are finite-dimensional because SU(2) 
(Spin(n)) is compact. " Deformative" generators of GL(3,M) (GL(n, M)) mix 
those representations and predict the existence of some excited "trajectories". 
And of course all characteristics of energy levels are primary and implied by the 
group itself. When one deals with elementary level, they are not interpreted 
as following from the rotation and deformation of something in the physical 
space. The natural question appears as to the possibility of such an interpreta- 
tion. There were some attempts, among other ones those by Barut, R§,czka and 
coworkers, of interpreting elementary particles as quantized rotating or rotating- 
and-deforming objects [11 [6]. Spin quantum numbers would not be then some- 
thing primary, but rather some aspects of quantized motion in the configuration 
space of internal degrees of freedom. But how to explain then the half-integer 
spin? There is some natural hypothesis based on two-valued wave functions. 
Incidentally, the idea was suggested long ago by Pauli and Reiss [S31 [SB] who 
claimed that the demand of one-valuedness imposed on wave functions in stan- 
dard quantum mechanics might be perhaps weakened and sometimes under cer- 
tain conditions rejected (cf. also [37]). This is the case in quantum mechanics 
of rigid bodies or homogeneously deformable rigid bodies. The classical config- 
uration spaces S0(3,K), GL(3,M) (SO(n,R), GL(n,]R) for n > 3) are doubly 
connected and it is quite natural to admit wave functions defined not on them 
but rather on their covering groups [H [6l [53l [58] . This idea has to do with 
the projective representations of groups. The only point is that to maintain 
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statistical interpretation in S0(3,R) (SO(n,R)) or GL(3,R) (GL(n,R)), one 
must assume some kind of superselection. Namely, the kernels of projections 
from the coverings onto original groups are two-element subgroups. There are 
two natural subspaces of wave functions: ones which do not distinguish the el- 
ements of the kernel subgroup and ones differing in sign there. Functions from 
two different subspaces cannot be superposed. But it is not clear if statistical 
interpretation must hold in original groups or in their coverings; in the latter 
case there is no restriction for superpositions. If proceeding along such lines, 
one obtains half-integer values of the angular momentum of quantized internal 
motion. It is interesting that, in a sense, something like this appeared also in 
our model of quantum torsional oscillator, even without working in the covering 
groups [68] . 

1 Degrees of freedom, kinematical and canoni- 
cal variables, symmetries 

We shall consider a system of material points, let us say physically a molecule or 
perhaps even the atomic nuclei, subject to translational motion, rigid rotations 
and homogeneous deformations. By this we mean that the current positions of 
constituent particles in the physical space, analytically speaking its Euler coor- 
dinates, are related to the reference Lagrange coordinates by the following 
formula: 

= X* + (^Va^^, (1) 

where a;* are coordinates of the spatial position of a distinguished point of the 
body, usually the centre of mass, and ip^ k refer to the relative/internal motion. 
Generalized coordinates are given by the system (a;% '■p^ k)] just these quantities 
are considered here as the functions of time. Generalized velocities are given 
then by the time derivatives i*, p^ Ki i-e., respectively the translational velocity 
and the system of internal ones. Sometimes it is convenient to use the quantities 

n j = LP 3, \l B ^ f t^P B ^ (fi i'^ijrB, {^) 

so-called affine velocities respectively in spatial and co-moving representations. 
If gij , rjAB are metric coefficients respectively in the physical and material spaces 
(Euler and Lagrange metrics) , then the corresponding doubled skew-symmetric 
parts, 

= n',-g'''gjin'k, n^B -^b^ ^n^B -v^^vBon'^c, (3) 

may be interpreted as angular velocities. Obviously, the upper case indices in g 
and T] refer to contravariant inverses of gij and tjab, i-e., 

g^'gk.^s^,, v^''vcB = 5^B. (4) 

Usually, although not always, one uses orthonormal coordinates in which 

ffifc =* Stk, VAB Sab- (5) 
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Remark: unlike in the skew-symmetric parts ^ in general are not (p- 
related, 

n', - n,' ^ p\ {^"^B - f^B-^) ^-^""j- (6) 

The equality holds if and only if ip is an isometry, i.e., 

Vab = gij^^A^^B- (7) 

Then f2*j, are respectively g- and 77-skew-symmetric and coincide with the 
angular velocity in spatial and co- moving representations. 

An important point is that the quantities fl^j, fl^B are non-holonomic veloc- 
ities or quassivelocities in the sense of Boltzmann. By this we mean that there 
are no generalized coordinates x^j or y^B such that Q^j, H^b would be their 
time derivatives. Quasivelocities are linear functions of generalized velocities, 
however, with configuration-dependent coefficients. 

In certain formulae it is convenient to use the co-moving components of the 
translational velocity, 

d^ = ^-'^y = ^-'^^. (8) 



The Green and Cauchy deformation tensors are respectively given by [221 [24] 

Gab = gioV^ aV^ B, Ci-j ^ rjAB'^^^^i'^"^^ j, (9) 

and their inverses are as follows: 

G = If lip jg\ C ^ ip Ar bV ■ (10) 

Warning: to avoid mistakes it is better not to omit the inverse label at C~^, 
G~^. Otherwise some confusions would be possible with the g- and 77-shifts of 
indices. And again it is typical that those are different things and except some 
special situations the following inequalities hold: 

C-H. ^ g^^^g^^Ckl, G-'^^ ^ V^'^V^'^GcD- (H) 

Having at disposal two pairs of twice covariant tensors {g,C), {ri,G), we can 
construct two mixed tensors 

= g^f^Ck, , G^B = V^^GcB (12) 
and the corresponding basic deformation invariants [22l [24] 

= Tr ) - Tr (c-^p) , p = 1, . . . , n. (13) 

The quantity Ip for any other value of p is a function of the above ones. This 
is a consequence of the Cayley-Hamilton theorem. For non-deformed config- 
urations the quantities G, G coincide respectively with r/, g. Sometimes it is 
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more convenient to use measures of deformations vanishing at the non-deformed 
configurations, so-called Lagrange and Euler deformation tensors p2[ 124) : 

E=\{G-fi), e^\{g-C). (14) 

Constructing from them the mixed tensors E, e, i.e., 

E^B^ri^^'EcB, =5''efc„ (15) 

we can obtain new versions of basic deformation invariants 

Tv(epY Tr(e*'), p=l,...,n. (16) 



They have the advantage of vanishing when there is no deformation; but of 
course there is nothing essentially new in them. They are some functions of 

®- . 

Transformation groups and symmetries are very important for our analysis. 
Linear spatial and material transformations act on internal degrees of freedom 
as follows: 

^^{Lao Rb) iv) = {Rb o La) [v) = A^B, (17) 

analytically 

[^V] ^ [A'.ip^LB'^K] . (18) 

If A, B are isometrics, i.e., 

5ijA'„i^-'n = gmn-, VKlB^ mB^ N = VMN, (19) 

then obviously Ip are preserved by (I17p : 

Ip{A^B)=Ipiip)- (20) 

this is just the reason why they are called deformation invariants. 

Affine velocities ft, suffer the following transformation rule under (|17p: 

n^AnA-\ h^B-^VLB. (21) 

Similarly 

[v^] ^ [A^,v'] , [v^] ^ [B-^"mv^'] ■ (22) 

Euler velocity field in the physical space may be expressed in the following way 
through the quantity fi^j, 

v\y)=v' + n'j{y^ ^x'), (23) 

i.e., is the gradient of y — ^ v{y) in affine motion. Similarly, ip^ a is related 
to the Lagrange velocity field. 
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Having in view analysis of equations of motion and first of all the quantiza- 
tion procedure, we must mention the basic concepts of Hamiltonian formalism. 
Canonical momenta conjugate to a;*, </?*a, or more precisely dual to generalized 
velocities i', <f^A, will be denoted respectively by pi, P^i- In certain formulae 
it is convenient to use the co-moving representation of translational motion: 

PA^PiV^'A- (25) 

Instead of P^i it is more convenient to use non-holonomic canonical momenta 
of internal motion conjugate to non-holonomic velocities 17'^, ^^b- We shall 
denote them respectively by S^s: 

S', - V'aP^,. S^B = P^^V'B = V-^^^^',^^B. (26) 

These quantities are referred to as components of afhne spin with respect to 
the space- and body-fixed axes respectively. They are Hamiltonian generators 
(momentum mappings) of transformations (fT7|) , (|2ip |^ |21 [Mj ■ Left and 
right acting (in the sense of (HJ)) rotation subgroups SO{n,g), SO {n,r]) are 
generated by the corresponding skew-symmetric parts, 

S', = E*,-S/' = S]^,-5,fe/S^, (27) 
V^B = ^^B-^B^ = ^^B-v''^'riBD^''c: (28) 

i.e., by the canonical spin S and vorticity V (in Dyson terms p[F). 
Warning (like in ©): V^b are NOT (!) co-moving components of S, i.e., 

V^B + ^-^^^S'j^'b. (29) 

The mentioned duality between non-holonomic canonical momenta and non- 
holonomic velocities is meant in the obvious sense 

E*,ri^; + p,v' = t^Bn'^A + PAV^ = P^^V'A + P^v' (30) 

for all virtual velocities. 

The non-holonomic character of j or T,^ b consists in that their Poisson 
brackets do not vanish; cf. below the formulae ((60)) . (|6ip. 

In practical calculations, especially ones concerning problems with high dy- 
namical symmetries, one uses special systems of generalized coordinates based 
on the polar and two-polar (singular value) decompositions. They are motivated 
physically by the classification of degrees of freedom and distinction between 
various modes of motion. 

Roughly speaking, the polar decomposition consists in representing tp as 
a product of linear isometry (orthogonal matrix) and positively definite sym- 
metric matrix in that or opposite ordering; so, there are two versions of this 
decomposition. 

More precisely, having two pairs of "metrics" (G'[(/7], ry), (C[(^], 5) for any in- 
ternal configuration Lp, we can express G[p\, C[(f \ through normalised orthogonal 
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principal axes (in the generic case those axes are unique): 

GM = ^A,F«M®F« (31) 

a 
a 

where -F'a[<p], fa[^] are vectors of orthonormal bases (respectively in rj- and g- 
sense) diagonalizing the deformation tensors G[ip], C[ip], and F'^[ip], /"[y] are 
elements of the corresponding dual bases. 

Let us mention that in this way the <y9-dcgrccs of freedom split into three 
subsystems: two fictitious rigid bodies (given by principal axes of the Green and 
Cauchy deformation tensors and deformation invariants, i.e., pure stretchings 
Aq). There is exactly one isomctry U[(p] such that 

U[^]Fa[ip] = faM, a=l,...,n. (33) 

It is clear that 

gijU[ipYKU[ipYL = VKL (34) 

and 

ip = U[ip]A[ip], (35) 

where A[ip] is ry-symmetric, 

mMA[^]^L = llLMA[ip]^K, (36) 

and positively definite. 

As mentioned, one can also write 

V = B[^]U[^], B[^]=U[^]A[^]U[^]-\ (37) 

where -B[<^] is ^-symmetric, 

gikB[ip]^j = gjkB[ip\\, (38) 

and positively definite. It is well known that the polar decomposition in both 
versions is unique. Analytically this is just the factorization into the product 
of orthogonal and symmetric-positive matrices. Diagonalizing the symmetric 
part one obtains the singular value decomposition, i.e., representation of y as a 
product of orthogonal, positive- diagonal and again orthogonal matrices, 

ip = LDR-^, i,i?e SO(n,M), £> = Diag (Q\ . . . , Q") . (39) 

Using more geometric language, i.e., making a systematic distinction between 

the material space (Lagrange variables), physical space (Euler variables) and 
M" (deformation invariants), we would write that 

if' A = L\D\R-^''a, (40) 
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where Z) : M" — »• M" is given by 

D\ = if aj^b, D% = Q''. (41) 

In some formulae it is convenient to use variables g", where 

Q'^ = exp {q-) . (42) 

The fictitious rigid bodies L, R, i.e., attitudes of the main axes of inertia have 
their own angular velocities and canonical spin variables. "Spatial" angular 
velocities are respectively given by 

Much more convenient are their "co- moving" representants: 

The corresponding canonical spin variables are respectively denoted by 

Q j, L, e b, T b- (45) 

Canonical momenta conjugate to Q", are denoted respectively by 

Pa, Pa- (46) 

Let us observe that Q'*, q", a = 1, . . . , n, are alternative measures of deforma- 
tion invariants. The two-polar splitting is not unique. The invariants or 
are indistinguishable "particles" on R and may be permuted; those permuta- 
tions must be accompanied by multiplying the matrices L, R on the right by 
appropriate orthogonal matrices having in each row and column only zeros but 
once ±1. More precisely, this finite non-uniqueness is a characteristic feature of 
the situation when the spectrum of D is non-degenerate, i.e., when all Q^-s are 
pairwise distinct, Q* if i ^ j, i.e., 

n (Q' - Q') ^ 0- (47) 

If degeneracy occurs, i.e., when the above product expression vanishes, then 
some additional continuous non-uniqueness of the right-multiplying gauge or- 
thogonal matrices also appears. An extreme special case is when all D*-s are 
equal, i.e., when the diagonal matrix D is proportional to identity matrix /„. 
Then L, R separately are not well defined and it is only LR~^ that is meaning- 
ful. The set of orthogonal right-multiplying gauge matrices is identical with the 
total orthogonal group. And there are of course intermediate situations between 
the non-degeneracy and the total degeneracy; all they are characterized by some 
continuous parameters. 
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It is clear that g, t coincide with spin and minus vorticity, 

q'j = T^L = -V'l. (48) 

The pairing between canonical variables and velocities (holonomic and non- 
holononiic) is given as follows: 

^p^q- + lTTigx) + ^TT{Td) = par + ^Tvigx) + ^TT(Tdy (49) 
Obviously, 

Par - PaQ". (50) 

Remark: do not confuse the above with translational momenta conjugate to 
i". When there will be a danger of confusion, from now on canonical momenta 
conjugate to i° will be denoted by p(tr)a. 

Deformation tensors are expressed by the following formulae: 

Gkl - VRsA'^kA^l, = gfe,B-i^',B-i',, (51) 

G^L = R"" aD\D'kR-^\, = L\D-^\D-^\L-'^' (52) 

In purely analytical matrix terms, when orthonormal axes are used and matrices 
of 77, g coincide with the Kronecker symbol, we can simply write 

G = = RD^R^ = RD^R-^, C = R-"^ = LD-'^L^ = LD-'^L-'^ . (53) 

One should mention what is the physical meaning of the above decompositions. 
Roughly speaking, U[ip\ in (l35|) . ([37]) describes rotational degrees of freedom 
and A[ip\, B[ip\ are two alternative descriptions of homogeneously-deformative 
modes of motion (equivalent respectively to the Green and Cauchy deformation 
tensors). But there are two important informations in deformation tensors: 

1. Purely scalar information about the stretching; in 71-dimensional space 
there are n independent ones (physically n = 3, sometimes 2 or 1). 

2. Information about how this stretching is oriented with respect to the 
space- or body-fixed reference frames. 

The first information is analytically described by deformation invariants, thus, 
finally by the diagonal matrix D. The "directional" information is equivalent 
to the knowledge of principal axes of deformation tensors G and C meant in the 
sense of metric tensors rj and i.e., normalised eigenvectors of G and G . And 
this information is contained in objects L and R. 

When some origin of the Cartesian coordinate system is fixed in the physical 
space, by convention the point y = 0, we can also introduce the orbital and 
total affine momenta ("hypermomenta" in terminology of Hehl and Ne'eman) 
[5U1 [5T1 [SH [75] with respect to that origin. We denote them respectively by 

A^- = x>(tr), , .r, = A', + S^- = xV(tr), + S^- . (54) 
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Their doubled skew-symmetric parts are respectively the "orbital" and "total" 
canonical angular momenta, 

L'j = a;>(tr)j - Xjp{try, J'^ = + S'j = x'p{tr)j - x^pitrf + S'j. (55) 

They are Hamiltonian generators of transformations 

'x^^VjX^, 'ip^^Vjip^K, (56) 

where L is orthogonal, 

gijL^kL^m — gkm- (57) 

Similarly, ([M]) are Hamiltonian generators for the non-restricted transformations 
with the general L, without the orthogonal constraints (|57|) . However, one 
must always remember that unlike the absolutely defined transformations (jl7p , 
(|56p are always related to some fixed origin in the physical space and so are the 
quantities ([55| . 

To finish with this kincmatical introduction, let us quote the basic Poisson- 
bracket relations finally following from 

{x\pitr),} = 6^,, W\,P''j} - S^aS^j, (58) 
{a;\P^,}=0, {x\ip^A}^0, {^'^,p(tr),} =0. (59) 

The most important of them read 

= J^S^ - (J^S^, [^^B,p{tr)c} = S^cp{tT)B, (60) 

[^''b,^^'d}^S^'b^^d-6^d^^b, {s'„E^b}=0, (61) 
{j^-,p(tr)fe} - {A^,p(tr)fe} = S^pitr),. (62) 

If some function F depends only on the configuration variables tp, but not on 
the generalized momenta, then 

{^,AM=.'^, ^''^ 

Non-quoted Poisson brackets do vanish. 

The system of Poisson brackets is very helpful when deriving equations of 
motion in the balance Hamiltonian form, 

f = WH). (64) 

where F runs over some complete family of Jacobi-independent functions on the 
phase space. 
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2 Inertial properties and canonical symmetries 



The summation of kinetic energies of constituents of our "molecule" results in 
the foUowing expression: 

T = T,. + Ti„, = -g, — — + -5,— , (65) 

where m is the total mass of the body and J^^ are, roughly speaking, the 
CO- moving components of inertial tensor or, more precisely, the second-order 
multipoles of the constant Lagrange distribution of matter in the space of La- 
grange coordinates, 

m = y dfiia), J^^ = j a^a^dnia) = J^^. (66) 

Here /i is the Lagrange co- moving distribution of mass, automatically constant 
in time; therefore, also m and J^^ are constant in time. 

Obviously, the formula ()65|) is valid under the assumption that are current 
coordinates of the centre of mass, i.e., that the origin of Lagrange coordinates, 
a^^ = 0, is chosen in such a way that the dipole moment of fi vanishes, 

= j a'U^i{a) = 0. (67) 

Otherwise in the kinetic energy expression T some interference, crossing terms 
bilinear in [dx'^ / dt,dip^ a/ di) would appear. For afFine bodies the vanishing of 
is equivalent to the equation 

y (y^-a;^)d/i(,,^)(y) = 0, (68) 

where pi(x,ip) is the measure obtained from /i as its transport by the mapping 
([1]). For non-affine bodies, when ([T]) is replaced by a general function (a^), 
things are more complicated, but there is no place here for their discussion. 
In any case, the condition = is not then equivalent to taking as a;* the 
coordinates of the current position of the centre of mass in the physical space. 

Incidentally, later on (the reasons will become clear after some additional 
remarks) it is convenient to express (|65p in the following form: 

r = rtr + rint = '^Gabv^v'' + ]^Gkl^'^ a^'^bJ^'' 

= "^g^jv'v^ + \9^,n\n^iJ[^f\ (69) 

where J['p]'-' are spatial, therefore time-dependent, components of the inertial 
tensor, i.e., 

JM'^' ^/'Vs^^''- (70) 
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Incidentally, (|70p may be written in the following form: 



J[vr = j (y' - x") {v' - x') dfi^^^^y (71) 

This kind of expression is applicable to general bodies and in the case of affine 
motion it becomes ([TO)) . Although J[ip] are components with respect to the 
incrtial frame, apparently more "touchable" than the co- moving (material) one, 
they are not very useful in dynamical equations. Let us remind that even in 
elementary problems of rigid body mechanics, we simply must use constant 
CO- moving components J^^ . Those are either calculated on the basis of our 
knowledge of the mass distribution in the body or somehow postulated, esti- 
mated a priori and testified on the basis of comparison of calculated results 
with experimental data. 

It is clear that (j65p is invariant under isometries acting in the physical space 
(parameterized by Euler variables) and under the group O {U, J) preserving the 
fixed inertial tensor J. If J is isotropic in the sense J^^ = Irj^^ , then we are 
dealing with the double orthogonal symmetry. 

If we use the above-mentioned polar decomposition (|35p and the correspond- 
ing "co-moving" angular velocity 

then the expression for the internal kinetic energy becomes 

Tint = -^Tr {AJ.AQ') + Tr (m^^) + ^Tr (^J, (^^^ j , (73) 

where Jrj is obtained from J by the 7y-lowering of the second index, 

J/b J^'^VCB- (74) 

Obviously, if orthonormal coordinates are used, there is no numerical distinction 
between matrices of J,, and J. The first term in (j73p is centrifugal one, the 
second one represents the Coriolis inertial forces, and the third one describes 
the dynamics of pure deformations. 

If we use the two-polar decomposition ([39]) and the body is doubly isotropic. 



jAB^j^AB^ (75) 

i.e., in orthonormal coordinates 

then the internal kinetic energy ([73]) can be rewritten as follows: 



Tint = -^Tr {D'x') - ^Tr [o'd'^) + /Tr [oxD^) + ^Tr [ ] , (77) 
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where the first two terms are centrifugal ones, the third one describes the Coriolis 
forces, and the last one is the kinetic energy of pure stretchings. If we perform 
the Legendre transformation for the Lagrangian 

L = r-l/(x,(p), (78) 

where T is given by (|77l) . and use the variables 

M:=-Q-?, N:=Q-T, (79) 

then ((77)1 becomes as follows: 

Although everything is in principle clear from the context, it is perhaps 
convenient and instructive to quote some explicit formulas and statements. 

Let L be a Lagrangian of some classical system with generalized coordi- 
nates ^'^ and generalized velocities Geometrically speaking, are 

adapted coordinates in the tangent bundle TQ, where Q denotes the configu- 
ration space, a manifold of generalized "positions". Legendre transformation C 
leads from generalized velocities to their conjugate momenta tt^; in geometric 
terms, it maps TQ into the cotangent bundle T*Q parameterized by (^'',7r^). 
Analytically it is given by 



(r,e-)-(e^-.)=(e',||). (81) 



dL 

In "usual" mechanics this mapping is invertible, more precisely, it is a diffeo- 
morphism of TQ onto T*Q. Inverting it we express generalized velocities in 
terms of canonical momenta; substituting this expression to the energy function 



^(^M^^A.) ^^M^_7,(^A.^^p) ^ (82) 



one obtains Hamiltonian as a function H : T*Q M, analytically, i/(^,7r). 
Later on one can proceed in Hamiltonian terms, following the formula (|64p . If 
Lagrangian has the potential structure ([75)) (sometimes called "natural" one), 
then obviously in ((8T|) L may be replaced by T itself. 

Some comments are necessary concerning the formulas (j79p . (|80[) . the more 
so that " non-holonomic" quantities X: "^j Qi ^ are used there. Non-holonomic ve- 
locities, or quasivelocities, introduced long ago by Boltzmann and many others, 
are given by expressions: 

LU^^LU^^iOt, (83) 

where depend on generalized coordinates ^ in an irreducible way, i.e., the 

Pfaff forms lj^ ^{^)d£_'-'' fail to be differentials. In many problems, appropriately 
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chosen are more convenient than typical examples are angular velocity (in 
space dimension higher than two) and our affine velocity. There is a conjugate 
concept of non-holonomic momenta, or quasimomenta, aA given by 

(TA^p^,a^'A, (84) 

where the uj- and cr-matrices are mutually reciprocal, 

uj^^af^B^S^B, a^Ai^^.^S^., aAio^ = 7Tj^. (85) 

The "non-holonomic" character of as is that their Poisson brackets do not 
vanish. Typical examples are angular momentum of rotational motion or our 
afhne momentum ("hyperspin"). And in general, when Q is a Lie group, then 
it is typical that uj, a are elements of its Lie algebra and co-algebra (dual of a 
Lie algebra). 

When non-holonomic quantities uj, a are used, Legendre transformation (|82p 
may be alternatively expressed as follows: 

(^^..^)-(e^a^) = (e^||), (86) 

where, obviously, L is to be treated as a function of ^, uj. The energy function 
E becomes 

E{tu;)=uj^^~L{^,uj). (87) 

Inverting Legendre transformation, i.e., expressing uj through a, one obtains 
Hamiltonian H (^,a). 

Just this procedure is meant in formulas ([77)1 - ([50]) . Configurations ^ are 
there given by the triples (Q, L, R) "parameterizing" the object ip. Non-holono- 
mic velocities we use are arrays {Q"' ,X°'b,'&°'bj , the kinetic energy is expressed 
by them in fTT)) : this is just our specification of the kinetic part of ((87|) . and our 
non-holonomic conjugate momenta are arrays (Pa, ^6, ?"&)■ The kinetic part 
of the Hamiltonian is given by (|80p , the total Hamiltonian contains in addition 
some potential energy term V [x; ip) =V {x; Q, L, R). Legendre transformation 
corresponding to (|77p is analytically given by the following expression of Pa, 
g^b, t\ through Q^ 

(88) 

I (^D^X + XD^ -2DdDy (89) 

/ (^D^d + - 2DxDj . (90) 

Let us remind that the matrix D is diagonal and Q"^ are its diagonal elements 
— deformation invariants. The last two formulas are a bit "symbolic" because 
X, 1? are skew-symmetric and their matrix elements are not independent. We 



P, 



Q b 



T b 



dx'a' 



9^ Wb] = 



[rb\ = 
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mean the following: In expression ([77| we replace Q"^, Xi by X + ^Xj 

+ W and calculate the corresponding variation of Tint, i-e., ^Tjnt up to terms 
linear in SQ"', Sx, Si!). After some easy calculations we obtain that 

<5Ti„t = PaSQ'' + ixr (gSx) + ^Tr (rSd) , (91) 

where the skew-symmetric matrices g, r are given just in ([55)1 - 1(511)) . Compare 
this also with ([50)1 . The space of skew-symmetric matrices is isomorphic 

with its own dual in the sense of the above trace expressions, and in this sense the 
last two formulas in ([88l) - ((90)) are meant, because the differential of a function is 
by definition given by the main part of its variation, i.e., one linear in variation 
of independent variables. The factors 1/2 in (|49p. (|88p - l|5D|) are due to the fact 
that because of the skew-symmetry of x, t every term in the traces in (|9ip 
occurs twice. In the physical dimension n = 3 (but neither in planar problems 
n = 2, nor in the "academic" situation n > 3), the above formulas may be 
written in a simpler and more familiar form, when x, Q, t are identified with 
3-dimensional pseudovectors (axial vectors); do not confuse them with the usual 
("non-pseudo") vectors: 

x\ = ~e\,r, ^% = -e^c^^ (92) 

Tb = e\^g,, ?\ = £%^re. (93) 

Obviously, e denotes here the totally antisymmetric Ricci symbol and the shift 
of indices is meant in the trivial sense of the Kronecker-delta metric. Then 
- (Plj) become respectively as follows: 



IQ\ (94) 



= (95) 

Ta = (96) 

STint = PaSQ^' + gaSr+TaSd-. (97) 



Solving (188l)-(|90l) with respect to D, x, one obtains ^ with M, N expres- 
sions defined in ([79]) . 

As seen from the very form of Tn- in (|65p . the Lagrangian given by Ttr alone, 
without the potential term V{x,ip), is non-physical in elasticity and condensed 
matter theory, because it predicts an unbounded expansion of the body. This is 
also seen in (j80p where the interaction between deformation invariants is purely 
repulsive. 

And at the same time there is some "aesthetic" drawback of (pS]) and the 
mentioned resulting formulae. Namely, in spite of the affine symmetry of degrees 
of freedom, the kinetic energy is not affinely invariant in the sense of dynamics. 
Kinetic energies, i.e., the metric tensors on the configuration space, invariant 
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under the left or right affine translations are respectively given by the following 
expressions: 

r = rtr + Ti„t = '^VABV^v'' + ^c^B''D^''An''c, (98) 

T = Ttr + Ti^t = '^g^Jv'v^ + ^n'j^n^,n\, (99) 

where C'^b'^d, T^^j'^i are both constant and symmetric in bi-indices {^b,'^'d), 
There are no metrics which would be simultaneously left and right 
afhnely invariant and non-degenerate. The reason is that the affine group is not 
semisimple and contains the normal subgroup composed of translations. 

There exist however some models of the internal kinetic energy invariant 
simultaneously under the left and right action of the affine group. They are 
given by the linear combination of second-order Casimir invariants: 

T,, = ^Tr (n^) + I (Tr 17)^ = ^Tr + I (Tr h)" , (100) 

where A, B are constants. 

The B-term is a merely secondary correction. The main structure is con- 
trolled by the A-term. It is not positively definite and this property of the 
kinetic energy might seem embarrassing. It turns out however that at least 
for the incompressible body the negative contribution to Tint may encode the 
dynamics, without introducing any potential energy into Lagrangian or Hamil- 
tonian. Such geodetic highly-invariant models may be often solved explicitly, 
analytically, in terms of some well-known special functions on groups. In some 
situations they may be solved in terms of exponential functions on groups, at 
least to some extent. 

There are also other interesting models where the total kinetic energy, in- 
cluding translational one, is invariant under the left-acting isometry group and 
the right-acting total affine one. There are also models of opposite properties, 
i.e., invariant under the left-acting affine group and right-acting isometry one. 
The model affinely invariant on the right may be interpreted as a very drastic 
discretization of the Arnold description of ideal fluids [5] . In Arnold theory this 
was the Hamiltonian system on the infinite-dimensional group of all volume- 
preserving diffeomorphisms; in our model this is the finite-dimensional group of 
affine volume-preserving mappings. Such models may be useful when describing 
molecules and droplets of nuclear matter. It is interesting that models affinely 
invariant in the physical space may be also realistic and convenient in condensed 
matter theory. The point is that due to strong interactions and strong concen- 
tration of matter, molecules may be non-sensitive (in their kinetic energy terms) 
to the metric g of the physical space; instead, they may "feel", e.g., the Cauchy 
deformation tensor C as a metric object underlying the structure of the kinetic 
energy. Incidentally, such situations are faced with, e.g., in solid state physics. 
The kinetic energy of electron is then based on the so-called tensor of effective 
mass, not on the usual metric geometry. Moreover, the tensor of effective mass 
may have various exotic properties, e.g., it need not be positively definite. 
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If we insist on positive definiteness, then it is interesting that when some 
phenomenological constants are appropriately chosen, then the left- or right- 
invariant affine kinetic energies of the total (internal together with translational) 
motion may have the positive signature. And at the same time they have the 
structure admitting rigorous analytical solutions. 

The kinetic energy invariant under spatial isometries but not under the larger 
spatial group and simultaneously invariant under all material afhne transforma- 
tions has the following form: 

T = J9^lV'v' + ^-9^J^\nh9^' + ^Tr (f^^) + | (Tr . (101) 

The kinetic energy invariant under spatial affine group and under material 
isometries (but not under a larger group of material, i.e., Lagrangian trans- 
formations) has the shape as follows: 

T = ^r^ABV^v^ + ^VAB^^C^'^DV'''' + (^') + | (xr 1^)' , (102) 

where in both formulae above m, /, A, and B are inertial constants. 

Obviously, the two last terms in both formulae (|10ip . (|102p are respectively 
equal to each other and different symbols fi, 51 are used only for aesthetic 
reasons. 

Let us remind the formula ([55]) and notice that poip , (|102p may be written 
in the following suggestive forms: 

T = ^Gabv^v^ + ^-Gab^^c^'^dG-^^^ + ^Tr (f)^) + | (Tr n " 

T = "^C^.v'v' + ^-C,,n\n^iC-^''' + ^Tr (Q^) + | (Tr . 

If we use Lagrangians of traditional potential forms, i.e., 

L = T -V{x,ip), (103) 

then the Legendre transformation may be written in any of the following equiv- 
alent forms, cf. also formulas (|8T])-(I87I): 



u ^ 9L dT ^, dL dT 

Mtr). - ^ - ^ ^ ^ = ^'''^ 

There are also various mixed possibilities. 

Inverting these formulae, i.e., expressing velocities in terms of canonical mo- 
menta and substituting them to T, we obtain the kinetic terms T of Hamiltoni- 
ans. We can admit also some potentials V{x, (p) and consider the total Hamil- 
tonians. However, it turns out that if the body is incompressible, the kinetic 
term alone may encode the dynamics of elastic vibrations. This resembles the 
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Maupertuis principle where the dynamics (including some hidden version of the 
potential energy) may be encoded in the metric structure of the configuration 
space, i.e., in some geodetic model [3]. 

The mentioned model of incompressible body predicts both the bounded 
motion (below some energy threshold) and unbounded one. The only obstacle 
comes from compressibility, which must be stabilised by some potential term 
preventing the body from the dissociation (splitting) or contraction (collapse). 

So, it is convenient to separate the isochoric motion from the pure com- 
pressibility. In the case of affinely-invariant dynamical models this is achieved 
by introducing on the real line M the centre of mass of the logarithmic deforma- 
tion invariants, 

q^-{q^ + --- + q"), (106) 
n 

and the corresponding conjugate momentum, 

p = Pi + ---+Pn. (107) 

These quantities have to do with the uniform dilatations. This is rather ex- 
ceptional phenomenon and should be prevented by some potential V{q). As 
mentioned, other rotational and deformative modes of motion may be described 
in a satisfactory way by some geodetic dynamical models. Nevertheless, some 
potentials are also admissible. The most realistic of them are superpositions of 
binary and dilatational models: 

where the second term corresponds to the shear-like part of the motion. 
The next quantities to be used are Casimir invariants built of S, S, i.e., 

C{k) = Tr (E'') ^ Tr (e'^) , (109) 

and of their trace-less parts, i.e., 

CsL(„)(fc) = Tr(a^-) =Tr(a'=), (110) 

where 

a = S - - (Tr E) I„, = E - - (Tr fl) I„. (Ill) 

n n \ / 

We need also the squared magnitudes of spin and vorticity, 

\\S\f = -^TT{S'), \\V\f = -lTr{V'), (112) 

but ATTENTION: those two positively-definite quantities, i.e., and 
are not equal to each other when some deformation occurs! 
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The model of T afSnely invariant in space but metrically invariant in the 
system of particle (Lagrangian variables) has the following form: 

^aff-„.ctr ^ (^^2^ + ^ (Tr S)' + i-llV^f , (113) 

where 

a = I + A, (i=-'- '-^ ^, M= J • (114) 

For the model of T metrically invariant only in space but affinely invariant 
within the body we obtain that 

^mcf-aff ^ ^ ^ (Tr E)^ + i-||5f . (115) 

In other words we have respectively 



^aff-metr ^ _C^2) + -C{lf + -\\V\\\ (116) 

^r"-^^ = ^^(2) + ^C(l)2 + i-||5f. (117) 
Separating the dilatational motion and shear, we obtain respectively 



^aff— metr _ 1 /n\ 1 „2 

- 2(/ + A)^^^(")^^^+2n(/ + A + ni?)^ 
/ 



2(/2 - ^2) 



(118) 



/y-metr— afF _ 1 ^ /9^ j 2 

" 2(/ + A)^''^(">^ 2n(/ + A + nB)^ 

^ "5f. (119) 



2(/2 - ^2) 



These expressions become the doubly affine invariant model when we put / = 0, 
i.e., 

^Sf = :r7Ci2) - — TT^ T^Cilf. (120) 

2A ^ ' 2A{A + nB) ^ ' ^ ' 

And similarly, they may be obtained from p20p when we substitute A ^ I + A 
and introduce in addition the mentioned ||V^|p-, ||S'||2-terms. 

Obviously, the geodetic system with Lagrangian given by the internal kinetic 
term (|120p is explicitly solvable in terms of matrix exponents. The correspond- 
ing geodetics are given by 

ip{t) — exp (fit) ipo — ifo exp (^fltj , (121) 

where ipo, fl, are constants and 

n ip-'^Qip. (122) 
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As mentioned above, it is assumed here that the translational motion is ne- 
glected or dynamically independent on the internal/relative one. It is clear that 
ipo is the initial configuration at t — and f2, O are constant (thus, automati- 
cally initial) values of affine velocities this is the reason for notation used in 
(|122p . Those initial conditions are completely arbitrary. The above general solu- 
tion (jl2ip is a priori obvious. However, we are interested mainly in dynamically 
relevant quantities 5*, pi, M°'b, -/V°6. Their time dependence may be in principle 
extracted from the exponential formula (|12ip . although it is not an easy task; 
some implicit function and inverse function procedures are meant. But on the 
level of differential equations for the relevant quantities g*, pi, Af^b, N"-}, this 
would be hopeless. There is also some qualitative advantage from using the ex- 
ponential description of geodetics. Namely, on its basis one can show that there 
are open subsets of bounded and unbounded trajectories ("below threshold" 
and "above threshold" motions) in the configuration space. And of course, it is 
so also on the level of the subfamily of state variables (q*, M°b, N°'b), because 
motion in L, i?- variables is always bounded; the submanifolds of corresponding 
degrees of freedom are compact. And the details of motion in L, i?-variables are 
not of large physical interest in problems like the scattering data, radiation and 
absorption. There are obvious quantum counterparts of these statements. 

For Lagrangians (|118p . (|119p the general solution is not any longer given 
by the exponential formula (|12ip . Only some special "stationary solutions" 
of such a form (obvious generalizations of stationary rotations in rigid body 
dynamics) do exist. In those solutions ipo is arbitrary, but f2, f2 must satisfy 
certain algebraic conditions. So, for (|118p . (|119p respectively the following must 
hold: 

= 0, [n,n<''^] = o, (123) 

where the rj- and transposes are respectively given by 

(n^^yB = v^'^'riBDn^c, {^'^Yj = g'^g.in'k. (124) 

Roughly speaking, the corresponding affine velocities are "normal" in the 77- 
sense for (|118p and in the g-sense for (|119p . 

Nevertheless, on the level of quantities (g\ M";,, A^°f,) the exponential for- 
mula is still useful and the time dependence of those variables may be, as for 
(jl20p . extracted from it. Again one must remember that in the exponential 
procedure based on (|120p . A must be replaced by {I + A). And everything 
concerning the bounded and non-bounded trajectories remains valid. 

Let us now quote the one-dimcnsional lattice aspects of the classical dy- 
namics [m [75] . They are interesting in themselves and are strongly related to 
the hyperbolic version of Sutherland integrable lattices (at least in the geodetic 
case). In the affine-affine version they have the following form: 

raff - 1 , 1 (M%)^ 1 iN%f 1 

int 2a4-^'^^32a^^sh2^ 32« ch^^ 2/3^ ' ^ 
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The reasoning which leads from (|100p to (|125p is a direct replication of that 
presented in formulas (|8ip - ([87l) . but now with the kinetic energy (|100[) used 
instead of ([77)) . 

In the explicit binary representation: 



Ana ^ ^ ' 32a ^ sh^^^^ 

32a fh 2 — 2- 2r7,ap 



or in a more apparent way: 



int 



W^'"^ 32A^sh2^ 



a^^b OT^b 2 

_ ^y-J^ + ^ ,2 ™ 

32A^ch^2l^ 2nM+nS)^" ^ ^ 

The lattice structure of the dynamics of deformation invariants is obvious. 
For the afRne-metric and metric-affine kinetic energies we have respectively 

^aff-metr ^ TTgf ^ / + A] + ^ ^ || , (128) 

^jnetr-aff ^ T^aff ^ J + ^] + __L_ || 5 1| 2 ^ (129) 

where 7-^ [A ^ / + A] denotes ([MH) with A replaced by / + A. 

It is interesting that the Casimir invariant C(2) has the following form for 
the incompressible body: 

CsL(„)(2)=Tr(a2) =Tr(a2) ^ i.^(p^_p,)2 (130) 



1 >p {M\f 1 ^ (iV%)' 
16 ^ sh^ 2l_£i 16 ^ ch^ 2l_£i ' 

aT^b 2 aT^b 2 

neither g nor p enters here (no dilatational contribution to dynamics). 

The mentioned possibility of encoding the dynamics of bounded elastic vi- 
brations in affinely-invariant kinetic energy forms (without using any potential 
energy) is explicitly visualized in formulas ()125|1 . (|126p . ()127|) . Namely, the neg- 
ative configuration describes the apparently exotic "centrifugal attraction" of 
deformation invariants. The total formulas describes the balance of attraction 
and repulsion which results in the existence of open families of bounded and 
escaping trajectories (open as subsets of the general solution, i.e., family of all 
trajectories). 
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It is seen that the deformation invariants g° behave hke indistinguishable 
"particles" on the real line R. On the quantum level they become some strange 
parastatistical " particles" . This follows from the very peculiar non- uniqueness 
of the two-polar decomposition, as described above. Deformation invariants, 
as seen in the above formulas, interact mutually and are coupled to each other 
via the quantities M"b, N^-f, which play the role of some "springs", respectively 
repulsive and attractive ones. These coupling "constants" are however variable; 
they are dynamical quantities which together with invariants g° satisfy the 
closed system of evolution equations obtained from ([55)) by substituting our 
Tint for the Hamiltonian H and our g", pa, M"-},, N°-b for the state variable F. 
Incidentally, the situation becomes particularly suggestive in the planar case 
n = 2, i.e., in the "Flatland" world as described by Abbott [1]. This has to do 
with the exceptional properties of GL(2,E) among all groups GL(n,K). These 
properties follow from the fact that S0(2, R), being a one-dimensional Lie group, 
is Abelian. Because of this, x = Xt'^ = '^^Q — Q — S,t^t^V and any 
of these matrices has only one independent component. Obviously, the same 
concerns the matrices M, N and we shall use the following symbols: 



(131) 



The two-polar decomposition takes on the following matrix form: 

LP = LDR-^ 



cos a — smck 
sin a cos a 





g2 



cos (3 sin (i 
- sin /3 cos /3 



(132) 



and as usual, in afhne models the exponential representation 

g'^ = exp (q'^) , a =1,2 



(133) 



is convenient. Canonical momenta conjugate to Q°- or g° are also denoted by 
the usual symbols Pa, Va- 

Canonical momenta conjugate to angular variables a, /3 are denoted by Pa, 
One can show that, obviously. 



1 -^1 c?" ql ^1 



£'^2 = = Pa, 



T^2 = T^2 = Pl3- 



(134) 
(135) 



It is clear that 



m^pp-Pa, n=pi3+pa. 



(136) 



It is also convenient to introduce the variables q, x, characterizing respectively 
dilatation and shear, 



1 / 1 2\ 

^{q +q ) , x = q 



The corresponding conjugate momenta p, Px are given by 



P=Pl +P2, 



^ {P2~Pl) 



(137) 



(138) 
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Then for (|125p we obtain that 



and similarly pi8|) . ()119|) become respectively 



^aff-mctr ^ 7:aff[^^J + ^^ + _j:^^ (140) 
^njetr-aff ^ -j:^ [A ^ I + A] + , (141) 

where, as previously, T^^f [A ^ / + A] denotes T^J^f with A replaced by / + ^. 
It is seen that the x-terms of the last three formulas may be written as 

'^nf[^] = ^+VZU^), (142) 
where the effective potential is given by 

y^ff X = 2 ^. 143 

™'"^ ' IGAsh^f le^ch^f ^ ' 

Obviously, a, (5 are cyclic variables, thus, p^, "pp are constants of motion and 
their fixed values m, n of (p^g — Pq), (p/3 +Pq) characterize the families of the 
X- motions in the general solutions. It is seen that if |m| < then at large 
distances the attractive negative term prevails; if |m| > then the repulsive 
term prevails for |x| — > 00. In the first case one deals with bounded vibrations 
in the cc- variable, i.e., in the mutual displacements of deformation invariants. 
As the first, repulsive term of V^^^x) is positively-infinite when x ^ 0, we 
are dealing formally with the characteristic shape known from the theory of 
intermolecular forces (cf., e.g.. Fig. 1 and Fig. 2 with x taken instead of K). 
If |m| > then l^'^„(a;) is purely repulsive. Thus, in the isochoric, SL(2,R)- 
ruled part of motion we are dealing with some thresholds |to| — \n\ separating 
the bounded (oscillatory) and non-bounded (decaying) motions. What concerns 
the (j-part, i.e., pure dilatations, motion is unbounded, just the uniform motion 
with constant velocity, thus, collapsing or infinitely expanding on the level of 
the variable exp((j). This is just the earlier mentioned necessity of correcting 
TTjjf by some extra potential V{ci) stabilizing dilatations (the potential well or 
something very quickly increasing with |g|) [19j. Here, in the exceptional case 
n = 2 everything is easily seen because M, N are now constants of motion. Ev- 
erything remains true for TT^^-motr g-metr-aff^ because Pa, are constants 
of motion. As well everything remains valid on the quantized level, where one 
deals respectively with the discrete and continuous spectrum of energy. 

All is so nice in affincly-invariant models. In traditional "d'Alembert" model 
([80)) only pure repulsion occurs in geodetic models and they are completely 
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non-viable without introducing of some potential term. In particular, in two 
dimensions rt = 2, (l80l) becomes 



Tint = ^ {P? + Pi) + {Q\ Q') , (144) 



where 



n2 



F^«'„ (Q\Q') = 2 ^ 2 . (145) 

"'"^ ^ 4/(gi-Q2)2 4/(gl+Q2)2 

Without the extra potential term, only the unbounded purely scattering motion 
of deformation invariants is possible. 

It is very interesting to deviate for a while from the hyperbolic Sutherland 
lattices (|130p to their trigonometric analogues. Namely, if we consider a dy- 
namical system on the unitary group U(n) and use again the analogue of the 
two-polar decomposition (|39p . 



ip^LDR-\ (146) 

where L,R £ SO(n,]R) are real-orthogonal and D is diagonal with complex 
unitary entries on the diagonal, 

D^a^expitq-), <7°eR, (147) 

then for the combination of Casimir invariants controlled by constants A, B we 
obtain the kinetic energy 



_ B , A„ /^.A B 

lint - 2 

where A > 0, B > and the Lie- algebraic elements 



^Tr {n') - I (Tr ^ -^Tr (n^) - f (Tr ll)' , (148) 



o = h = ip-^^ ^ ^p-'^nip (149) 

are skew-Hermitian. 

The corresponding kinetic energy becomes then in the Hamiltonian repre- 
sentation as follows: 

- 2A^^''^ 32^^sin2 2l^ 
Separating explicitly the dilatational ((7,p)-variables we obtain 

a^b a^b 2 

32A^^cos2^ 2n(A + nB)^- 
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The Casimir invariant C(2) for SU(n) takes on the foUowing form: 



a^b a=ib ^^"^ 2 a^b 2 

Let us notice an important role of the negative contribution to (|130p . It 
describes attraction and stabihses the vibrating regime of incompressible bodies. 
All contributions to (1152^ are positive. However, it does not matter in view of 
the circular topology of the g-variable and all g^-ones in general. Because then 
it is impossible to distinguish between repulsion and attraction. 

It is very interesting that in (fng|) the characteristic l|5'|p- and 

terms appear. Those invariants of spatial and material rotations are very sug- 
gestive and resemble certain formulae appearing in the description of Raman 
scattering, rotational and deformative excitation of nuclear matter, especially 
on the quantized level. 

Let us also mention, one can suspect some physical interpretation in collec- 
tive models, where both the ||5'||^- and ||y|p-terms appear. They would contain 
highly symmetric afhnely-invariant expressions controlled by A, B as above and 
in addition two orthogonal terms restricting the symmetry group to isometrics in 
the physical and material space (in Euler and Lagrange variables respectively) . 
In terms of velocities they would be given by the following phenomenological 
formulae: 

Tint = ^g^kg^'n^jn^ + ^VKLV^'''n''M^''N + ^Tr (o^) + | (Tr 
Obviously, the third and fourth terms may be as well written as 



(153) 



^Tr(^^2) + |(Tr^7)^ (154) 
The first term may be alternatively written in the following form: 

^-^GKLG-'^'^'n^Mn^N. (155) 

Similarly, the second one may be written as 

jakG^'^'n'^n^. (156) 

Every of these forms may appeal to some intuitions and may be suggestive when 
properly read out. 

It is clear that after the Legendre transformation based on (|153p we obtain 
the following kinetic term of the Hamiltonian: 

■lint — -;=gikg 2^ I + -^IJkLV 2^ M N 

III 2I2 
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where Ii, I2, A, B are some inertial constants built algebraically of Ii, I2, A, 
B. Let us also mention expressions similar to (|154[) . (|155[) . (|156[) . 

It is also clear that the corresponding expressions for (|157p will contain on 
equal footing both terms proportional to WSW^ and Let us also mention 

that the dynamics for g", pa, M°f,, N°'b is closed (at least when the potentials 
V do not exist or depend only on the deformations invariants q^, . . . ,q", first 
of all on q separated from g" — q^). Equations of motion written in terms of 
quantities Pa, N^-t are autonomous and independent of L, i?-variables. 

They have the following form based on the Poisson brackets: 

% ^ ^i' 

= {P..«) = -||. (159) 



dt 



where 

= {(z^Ar4 = {p„,i\/r4 = {g^7v^4 = {pa,iv-4 = o. (i62) 

To obtain the Poisson brackets for Af , N one must use the brackets for r 
from which M, A'^ are built. In turn, the brackets for g, r are implied by the 
commutation relations in the Lie algebra SO(n, R)' for SO(n, R) because p, r are 
Hamiltonian generators of transformations acting on L, R in (|40p . respectively, 
the following ones: 

[L^] 1^ [L\A\] , [i?^,] , (163) 

where A,B ^ SO(n, R). Those Hamiltonian generators are related to spin and 
vorticity variables as follows: 

= L-^\S\jL\, t\ = -R-^-KV^LR'-b- (164) 

Therefore, after the Kronecker-delta shift of indices we obtain that 

{Qab,Qcd} = dadScb - QcbSad + ddbSac - QacSdb, (165) 

{Tab^T'cd} = TadScb - Tcb^ad + Tdb^ac - Tac^db, (166) 

{dab,Tcd} = 0, (167) 

and finally, 

[Mab, Afc4 = [Nab, N^d} = McbSad - Madkc + MaJdb - MdbSac, (168) 

{Mab,N^d} = N,b6ad-NadSbc + NacSdb-NdbSac- (169) 
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Solving (in principle) equations of motion (|158p . (|159p . (|160p . (|16ip and per- 
forming partially the inverse Legendre transformation, one obtains that 

so, one can find (in principle) the time dependence of x, •& and then the time 
dependence of L, R may be found by solving the following equation: 

dL ^ clR ^ 

3 Quantum models 

We are interested in studying physical phenomena in the micro- and nano-scale, 
where the quantized theory must be used. 

Roughly speaking, for quantum models the configuration space of inter- 
nal/relative degrees of freedom may be identified with (GL(n, M)), a bit more 
rigorously with (LI([/, V)), where LI([/, V) denotes the manifold of linear iso- 
morphisms of U onto V (usually the volume-preserving ones). Kinetic energies 
are usually based on some underlying metrics. 



TMdy'-^dy^ T=-r^.(y)^^. (172) 



The corresponding Hilbert space of wave functions (Q^fJ-) consists of square- 
integrable functions on the configuration space Q endowed with the canonical 
measure /i built of the metric tensor T. The corresponding scalar product is 
given by 

(«'i,*2) = / *i(y)«'2(y)dAi(2/), (173) 



where 

dti{y) = ^\det [r^,(y)]|V---dy^, iV-dimQ. (174) 

Fortunately, our metrics have certain invariance properties and because of this 
the measures fj, are what mathematicians call Haar measures [3^. Those are 
invariant under the left (Euler-spatial) and right (Lagrangian-material) trans- 
lations, cf. (fT7|) . If we take into account the translational degrees of freedom, 
then the corresponding measure a is given by 

da{x,(p) — (det iy9) " ^ dx^ ■ ■ ■ dx"d(p^ i ■ ■ ■ dip"n 

= {detip)-^dX{ip)dx^ ■■■dx'', (175) 
dX{ip) = {detipy"dip\---dip'\,. (176) 

The two-polar (singular value) decomposition implies that 

dXiip) = dX{L, g, i?) = J||sh {q' - q^) \ dfi{L)dfi{R)dq^ ■ ■ ■ dq"-, (177) 
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where fi denotes the Haar measures on the compact, connected and simply 
connected orthogonal groups SO(n,R). 

If we wish to consider the incompressible bodies, then the dilatation factor 
must be cancelled by the corresponding Dirac distribution, 

dXsLin){v) = n l^li - 9') I dfi{L)dfi{R)S (qi + • • ■ + g") dq' ■ ■ ■ dq^. (178) 

When quantizing the "d'Alembert-like" models ([50]) . it is more convenient to 
deal with the usual Lebesgue measure € on the manifold of internal degrees of 
freedom, 

d£{ip) = dip\---dip"n. (179) 
Then, if we use the two-polar splitting, we have that 

de{L, Q, R) = Pe (Q\ . . . , Q") dfi{L)dn{R)dQ^ ■ ■ ■ dQ'^, (180) 

where 

Pi{Q\...,Q'')^l[\{Qn'~{Q'f\^ll\(Q' + Q')(Q''-Q')\- (181) 

If translational degrees of freedom are explicitly taken into account, then in 
analogy to (I175|) we have that 

da{x, if) = d£{(p)dx^ ■■■ dx" = di (i, Q, R) dx^ ■ ■ ■ dx". (182) 

It will be also convenient to write the Haar measure (|177p in the abbreviated 
form analogous to p80)) . p8T)) : 

dX{ip) = dX{L, q, R) = Px {q\ . . . , g") dfi{L)dfi{R)dq^ ■ ■ ■ dq" , (183) 

where now 

Px{q\...,q'')^ll\sh{q^-q=)\. (184) 

Obviously, the measure A is invariant under all transformations ([T7| : simi- 
larly, a is invariant under (|17p accompanied by all afline mappings acting on a;* , 
i.e., on translational degrees of freedom. Unlike this, the Lebesgue measures i, 
a are invariant only if the mentioned mappings are restricted to isometrics. 

In general, the procedure of Schrodinger quantization in a Riemannian man- 
ifold {Q,T) begins from introducing the operator of kinetic energy, proportional 
to the Laplace-Beltrami operator pTl I41j . 

where denotes the operator of the covariant differentiation in the F-Levi- 
Civita sense along the /i-th coordinate axis. 
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The operators {h/i)'V^ and — (fi,^/2)A are formally self-adjoint, i.e., satisfy 

[mill] 

(A^-il^a) = (*i|A*2), (186) 
if 'i'l, are confined to some dense subdomain of L^(Q,/i) consisting of suf- 
ficiently smooth functions. Being differential operators, they are evidently un- 
bounded. They are formally self-adjoint, because the Levi-Civita parallel trans- 
port does preserve the Riemann measure fj,. 

However, in general such a procedure would be extremely strenuous. It 
would be very difficult to avoid mistakes and even if avoiding them we would 
obtain some rather obscure, non-readable expressions. Fortunately, differential 
operators generating left and right regular translations in the configuration space 
enable one to simplify the procedure in a remarkable way. 

Namely, it may be easily shown that transformations of wave functions in- 
duced by the argumentwise action of (fT5)) are unitary in (LI([/, V), A), just 
due to the left and right invariance of the measure A. Similarly, the usual 
vector translations in the physical space, just as (fT8| themselves, are unitary 
in L'^{Q,a). On the other hand, there is a unitary failure in {Ll{U,V),£), 
L^((5,a), unless the argumentwise action of transformations on "if will be ac- 
companied by certain multiplicative correction factor built of the determinants 
of A, B in 

Because of all that, the operators 

= -^""K^n-, (187) 







^ b 






I 


B 





dip 



K 



h f) 

a " ~1A b -^a /ICON 

B = ---P B r. „ =y aP B^ b (188) 

t I Oip°- A 

are " formally Hermitian" . The " formally anti-Hermitian" first-order differential 
operators L°f,, R'^s are infinitesimal generators of the mentioned unitary groups 
in L^((5, a). Obviously, the operators of spin and vorticity, i.e., 

S% - S% - Sb"^ = S% - g''-5MS''c, (189) 

Y^B = S^B-SB^ = S^B-?y^Si3i?S^c, (190) 

or rather their (i/?i)-multiplies, are infinitesimal generators of "spatial" and 
"material" rotations. S"f,, S'^b are operators representing affine spin with re- 
spect to the space- and body-fixed axes respectively. An important point is that 
no problems of ordering operators appear here. Namely, just due to the geo- 
metric interpretation of operators as generators of well-defined transformation 
groups, the ordering is exactly like in (|187|) . (|188|) . 

Instead of fighting with formulae like (|185p we simply write the following 
well-defined expression for the model of internal kinetic energy affinely invariant 
both in the Euler and Lagrange variables: 



rpaff— aff ytz 

" 2A ' ' 2A{A + nB) ' ' 



1 „, B 
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just the automatic replacement of S'j, by the operators S'j, in the 

corresponding classical formulae. 

For the models only isometrically invariant in Euler variables and affinely 
invariant in Lagrange ones we obtain obviously the following operator of kinetic 
energy: 

TSf ^-'^'^ = y^ug''^'j^\ + + (192) 

11 2A Zti 

where 

7^^, 1=^, j,^Jl±mL±A±j}^, (193) 

Obviously, the kinetic energy operator aflinely invariant in Euler variables 
and only isometrically invariant in Lagrange ones can be obtained in the follow- 
ing dual form: 

^aff-metr ^ l^^^^^^CD^A^B ^ ^ is^fiS^^ + -Lj^^a'^^'b, (194) 

21 2A 2B 

with the same meaning of inertial constants (|193p . 

The accompanying expressions for the operators of translational kinetic en- 
ergy are given respectively by 

TtT'^^^^ = ^5''-'p(tr).p(tr), = ^G-i^^p(tr)^p(tr)B, (195) 

Ttf-"*^*^ = ^C-'''p(tr)Mir), = ^rj^^p{tv)Ap{tv)B, (196) 

where, let us remind, p(tr)i, p(tr)yi are formally Hermitian operators of linear 
(translational) momentum as expressed in spatial (laboratory)/material (co- 
moving) terms: 

On the classical level we used the logarithmic dilatational invariant q (|106p 
and its conjugate canonical momentum p (|107p . In quantized theory this mo- 
mentum is represented by the following formally Hermitian operator: 

P = ^f (198) 

It is interrelated to the shear parts of the afhne spin (deviator) through the 
following formulae: 

s%-S%--pJ%, s^B = S^B - -P<5^s, p = S^-S^A. (199) 

n n 
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Just like on the classical level, one can perform a partial separation of shear 
(incompressible motion) and dilatations. Expressions for the operators of inter- 
nal kinetic energy become then as follows: 

rpmctr— aff ^ (<r\\ 



1 



Taff — mctr ^ (<r\\ 

int - 0/r I /1\ *-"SL(n)l^J 



where on the quantum level we mean that 

CsL(„)(fc) =s%s^...s%s^ = s^BS^c---s«ss^A, (203) 
(k factors) and 

||S||2 :=-is%S''a, ||Vf :^-iv^BV^^. (204) 
Sometimes, however, it is more convenient to write simply that 



1 . 1 o 1 



1 1 o 1 



Trfr'"''*'' - ^<^(2) + -p2 + -||Vf, (207) 

where a, /3, are given by (|114p and again the operators C(fc) are built according 
to the Casimir prescription, 

C(fc) = S^S", • • • = S^bS^c • • • S^sS^A, (208) 

fc multiplicative factors. 

The quantized version of the model (|153p . p57p will be based on the following 
operator of kinetic energy: 

Tmetr— motr „j'v* V*-' _L ^ „ „-,l^IN^K 

Zii 21 2 

+ j=^^L^^K + ^^''K^\. (209) 



In all the above expressions for T there is no problem with the ordering of 

J' 



operators just due to the geometric interpretation of the operators S'^, S^s 
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as generators of some groups of unitary transformations. And it is easily seen 
that all the resulting expressions for T automatically are formally Hermitian. 

What concerns (|209p itself, this model is spatially and materially invariant 
only under the isometry groups. However, it has a nice structure because it is 
a superposition of two affinely-invariant terms and two additional ones which 
restrict this symmetry to a weaker one, namely isometric. Manipulations with 
phenomenological inertial constants enable one to control somehow those sym- 
metry properties. After some calculations performed on (|209p . we obtain that 

Tmctr-metr ^ 1^(2) + ^ p2 + 1 ||S||2 + ^ ||V|p, (210) 

where a, /?, Ji, v are some constants built of /i, I2, A, B, in analogy to (|114p . 
It may be also convenient to represent the above expression as 



^mcf-netr ^ lCsL(„)(2) + + ^||S|P + l||Vf , (211) 



where a, 5, c, d are some new constants. In any case, the formula (j21ip may 
be also postulated as something primary, without the intermediary step (j209p . 
just as a natural generalisation/unification of (|206p . (|207p . 

Just as in the classical theory, spin and vorticity operators may be expressed 
in terms of their components with respect to bases co-moving with the L- and 
i?-gyroscopes, 

r\ = L-'\LhS'j, t\ = -R-^"'KR''by'^L, (212) 

the ordering of variables as indicated. Due to the orthogonality of L, i? it is 
clear that the following holds for the "magnitudes": 

f%f^, = S',S^„ tV^-V^iV^K. (213) 

In geodetic cases, or with potentials V (g^, . . . , g") built of deformation invari- 
ants, S'j, 'V^ B are quantum constants of motion, i.e., they commute with the 
Hamiltonian 'H. — T + V . It is not the case with r, t, however their squared 
magnitudes, being equal to those of S, V, are also constants of motion. 

Just like in the classical theory, in certain quantum expressions it is conve- 
nient to use the following operators: 

M% - -v\ ~ t%, N% - r\ - t%, (214) 

which enable one to perform a " partial diagonalization" of the kinetic energy. 

It is clear that for all geodetic models or more general dynamical models 
with potentials depending only on deformation invariants, the eigenvalues of 
orthogonal Casimirs of spin and vorticity 

l|Sf - ||ff , ||Vf ^ ||tf (215) 
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are "good" quantum numbers. In the physical three-dimensional case they are 
given respectively by [83l [Ml [85] 

C{S,s) = h^s{s + l), C{V,j)=tej{j + l), (216) 

where s, j are non-negative integers, 

5 = 0,1,2,..., j = 0,l,2,... (s,je{0}UN). (217) 

It may be also shown |65l [7ll [72l [23| that s, j may be non-negative integers and 
half- integers, 

s = 0,i,l,|..., j = 0,i,l,|... (^s,je{0}uiN^, (218) 

with the condition, however, that s, j are simultaneously integer or half-integer, 
i.e., (j — s) is an integer: 

=0,1,2,... (Ii~s|e{0}UN). (219) 

This has to do with admitting some special kind of multivalued wave functions; 
the procedure suggested among others by Pauli and Reiss [S31 [5S] . The configu- 
ration space of internal degrees of freedom, originally identified with GL^(3,R), 
is then replaced by its universal covering manifold, i.e., the universal covering 
group GL+(3,R), which, by the way, is not a linear group (it does not admit 
any realization in terms of finite-dimensional matrices). 

There is an interesting message of formulae pOTI) . (|202l) . (|206)) . (|207)) . (|2T0l) . 
(I211|) . concerning the spectrum of radiation of objects described by Tint as 
Hamiltonians of internal motion. More generally, this applies also to Hamilto- 
nians of the following form: 

H = T + y(qi,...,g"), (220) 

i.e., ones with potential terms depending only on deformation invariants (scalar 
stretchings). Those spectral peculiarities appear, in particular, in phenomena 
like the Raman scattering, when the absorbed light gives rise to the internal 
excited states which decay through radiation proving the splitting of internal 
energy levels. This splitting is imposed onto some system of background lev- 
els corresponding to the spectra of the first two terms in ()20ip . (I202p . (I206|) . 
(|207p . (|210p . (|21ip and is described by the terms proportional to operators ||S|p, 
II V p. As we know, in the three-dimensional physical space these operators have 
spectra (|2T6l) . (|2T7|) . under certain conditions (|2T8l) . (|2T9ll . 

In expressions (j201|) . (|206p we easily recognise the rotational Raman split- 
ting controlled by the quantity ft^s(s -I- 1). These terms correspond exactly to 
excitation of rotations in the physical space. 

From this point of view the models (|202p . (|207p describe something else, 
although the splitting has again the structure h^j{j + 1). But this is not the 
quantized rotation. Instead, it is some part of the quantized dcformativc motion. 
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i.e., some aspect of quantized deformations. The "rotational" expression + 
1) is simply due to the rotation of squeezing, rotation of the deformation tensor. 
So, in spite of the + l)-structure this is not any rotation of the "molecule" 
in space, this is rather something like the rotation of some external factors 
suppressing the " molecule" . 

In (|210p . (|21ip one deals simultaneously with both aspects: the quantized 
?i^s(s + l)-rotation in space and the quantized h^j{j + l)-controlled deformation 
process. This might be something realistic, because in spectra of some micro- 
objects one observes splittings of the h'^j{j + l)-type which cannot be interpreted 
as a quantized rotation in space. 

In nuclear physics there appear terms of the type h^I{I + 1), where / is 
an isospin. It is so even in elementary particles, where the mass formula of 
eight-fold way, obtained by Gell-Mann and Okubo for hadrons reads: 

M^a + bY + c(^I{I+l)-W^^ , (221) 

where / denotes the isospin quantum number, Y is so-called hypercharge, and o, 
6, c are constants. This is particularly remarkable when we consider the model 
(|150p . (|15ip . where GL(n, M) is replaced by U{n) is such a way that deformation 
invariants exp{q°-) are " compactified" to exp{iq°-). 

It seems that the invariance structure and symmetry groups are so fun- 
damental for dynamics that they may lead to quite similar models in rather 
mutually remote areas of physical phenomena. 

Incidentally, let us mention that the term linear in Y may have an analogue 
within our treatment, and namely if we admit in the formula for the kinetic 
energy some terms linear in generalized velocities. The only geometrically in- 
variant ones are those proportional to Tr = Tr i.e., proportional to the 
operator p on the quantum level. Apparently, the terms proportional to veloci- 
ties might seem exotic. Let us remind, however, that they appear in analytical 
mechanics of charged particles moving in the magnetic field. We did not con- 
sider above such models with kinetic terms linear in velocities, nevertheless, it 
may be easily done. 

Let us quote for the sake of completeness some formulae concerning the 
quantum description. In n dimensions our wave functions may be expanded in 
the series 

N(a) N(p) 

v&(^) = *(L,g,i?)= ^ ^ ^ P^„(L)/:f (g)Pf, (222) 

with the following meaning of symbols: 

• is the set of equivalence classes of unitary irreducible representations of 
SO(n,M), 

• N{a) is the dimension of the a-th representation class; as SO(n,R) is 
compact, N{a) is finite. 
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This follows from the Peter- Weyl theorem [HI [Si] applied to SO(n,K). 

In the physical case n = 3, is in principle the set of non-negative integers, 
a, (3 are, just as above, denoted by s,j = 0, 1, 2, . . ., N{s) = 2s+l, N{j) = 2j+l, 
and the indices (m, n), (fc, I) are considered as jumping by 1 from —s to s and 
from —j to j respectively, and P^, T)^ are standard expressions for unitary 
irreducible representations of S0(3,]R). As mentioned, according to certain 
ideas by Pauli, it is possible to admit some two-valued wave functions or 
more rigorously, wave functions defined on the covering group GL(?i,K). This 
group is nonlinear, i.e., non-realizable in terms of finite matrices. 

As mentioned, a, f3 (s, j) are "good" quantum numbers, so it is often con- 
venient to use just the reduced amplitudes 

N(a) JV(/3) 

^^I'liv) = "^"JiiL, 9, i?) = E E K^n{L)f:^{q)K {R-') ■ (223) 

n=l k=l 

In the physical case n — 3 this becomes 

^Zi^) = ^ZiL,q,R) -EE T^ZAL)a{q)Vi, . (224) 

n= — s k=—j 

As mentioned, when two-valued wave functions are admitted, then S0(3,R) is 
to be replaced by SU(2) and in the above series only such terms appear that 
(j — s) is an integer. Obviously, the following eigenequations hold: 

llSf*^, = \\vf^Z = ^'^{^ + mZ. (225) 

m'^z = \\tr^z-f^'jiJ + ^)'i'Z- (226) 

Traditionally one uses the convention that m, / are related to the eigenvalues of 
S3, V3, the third components of spin and vorticity: 

Sa*^, = hm^Z, Va*^, = M^Z- (227) 

Similarly, when the values n, k in the superposition (|224p are kept fixed and we 
retain only the corresponding single term, for the resulting reduced amplitudes 
we have that 

h'^Z = nn^Z , ts^Z = hk^Z ■ (228) 

In three dimensions, when are well-known functions found explicitly by 
Wigner, this means that the dependence of ^E" on "angular" variables L, R is in 
a sense explicitly known. And the action of differential operators occurring in 
our formulae is "algebraized". In the two-dimensional space, when n = 2, this 
is just the expansion into Fourier series. In three dimensions we have 

S^*"^ = S'^",,*"^ V^ij^-^^ = ^"^S'^^B, (229) 

where , are matrices of angular momenta indexed by s, j. In the academic 
general case s, j would have to be replaced by some labels a, f3. If n = 3, then 
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S", are standard matrices (2s + 1) x (2s + 1), (2j + 1) x (2j + 1) which are 
exphcitly known and quoted in any textbook on quantum mechanics, e.g., one 
by Landau and Lifshitz [3S] (cf. also [Ml [55] V 
In exphcitly matrix terms we can write that 

^'"^(L, q, R) = V°'{L)r'^{q)V^ {R'^) . (230) 
Differential operators f , t are algebraized as follows: 

r^^^^J°'f> = V°'{L)S°'''br'^{q)V^ {R'^) , (231) 
t^*"/^ = V°'{L)f"'^{q)S'^''bV'^ {R-^) . (232) 



In certain formulae it is convenient to use the symbols S'^'^b, S^°'b, where 



g»a^jal3 ._ s°"'bf"f^, S^'^bP'^ := P'^S'^'^b- (233) 

The Schrodinger equation reduces then to the infinite system of multicomponent 
Schrodinger equations for the reduced amplitudes given by the system of N{a) x 
N{P), i.e., physically (2s + 1) x (2j + 1), g-dependent matrices f"^{q). The 
scalar product (4'i|^'2) for the complete wave functions may be expressed in 
the following way through the one for reduced amplitudes: 

(*.i*^> - E^s^/Tr(/r» ..,")) 

X Pi{q' q")<lq^---iq", (234) 

where Pa (9^, ... , 9") is given by p84p and N{a), N{(3) are dimensions of irre- 
ducible representations a, /? £ 0. They are finite because SO(n, M) are compact 
groups. 

Similarly, for the d'Alembert-type models one uses the representation: 

(*ii*2) = E^^^^(^/Tr(/r'^(Q^...,Q'V2"''(Q^••■^^ 

X (Qi,...,Q")dQi---dg", (235) 

where Pf (Q^, . . . , Q") is given by pSTI) . 

The operator (|19ip for the doubly invariant kinetic energy may be written 
as follows: 



Taff — aff ' " T-k 

int " ^Wa^>' 



2 



2 A 2A{A + nB)dq 
^V^^-— (236) 

a,b '^^^ 2 a,b 2 



where the differential operator is given by 



^~P^^d^^d^^2^^diq''f^^ dq'' ^ ' 
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Let us observe that unlike to what might be naively expected, the operator 
Da involving differential operators d/dq"' is not the usual R" (physically B?) 
Laplace operator. One can reduce it to such a form by modifying the dependent 
variable, 

$ ^/p'x'^, (238) 
but the price is that an artificial amended potential appears: 

hi I ^ fdPxV , X 

^ = -2APF + 4ip;:?(a^) • (2^9) 

As mentioned, the time- independent Schrodinger equation, i.e., eigenequation 

H* = (240) 

splits into the infinite family of multicomponent amplitudes /"^ (q^, . . . , g") 
involving, however, only n independent variables q^, . . . , g" instead the primary 
r? ones ^p^ a, 

jjQ/3 jap _ j^ap jaP (241) 

where the reduced Hamiltonians H"'^ have the following form: 

■trap rap _ D . f"^ 

J 2A 

2 ,< ,x 2 



V / ^ \ V / fa0 



a,b 2 a,b "^^^ 2 



r 



fc2 rj q2 /"q/S 

^ +F(g\...,g")r^ (242) 



2A(A + nB) 



where the meaning of symbols 5'"";,, S^°'i, is like in (|233p . The potential energy 
is assumed to depend merely on logarithmic deformation invariants q^,. . . ,q". 
As mentioned, on the level of incompressible motion, i.e., for SL(n,R), it is 
possible to remain on the purely geodetic level. Then it is sufficient to admit 
potentials V{q) depending only on the dilatational parameter q (|106[) . The prob- 
lem splits then into the geodetic one on SL(n, E) and one-dimensional quantized 
oscillations in the g-variable. As usually one deals with almost incompressible 
(almost isochoric) motion, as for V{q) some simple phenomenological model 
may be used, e.g., some potential well or "steep" oscillator. In metric-affine and 
affine-metric models, (|242p is respectively modified by the following terms: 

-^Cia,2), "^C(/3,2), (243) 

based on second-order Casimir invariants for the orthogonal group SO(n,R). In 
the physical case n = 3, they become [52] 

^s(s + l), + (244) 
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with the previous meaning of symbols s, j. 

In more hypothetical metric- metric models (|210p . (|21ip . we are dealing with 
the following additive correction term: 

+ + + (245) 

Finally, let us quote a rather not very useful, although seemingly "more 
familiar", d'Alembert model: 



2/ {Q'^-Q'f 



1 ^ Cs^b + S^b)" 

YiT. \Q.^Qbf r' + v{Q\-,Q-)r^ (246) 



where 



^^-p,2^dQ^P^dQ-^-2^Jl^ + T,l.[dQ^) ■ (247) 

Let us note that for (|246p both the classical and quantum geodetic models would 
be completely non-physical and the use of some relatively general potential 
V [Q^, ■ ■ ■ , Q") is absolutely unavoidable. 

It is worth to mention that in planar problems, when n — 2, there exists 
a wide class of models which are both qualitatively physical and integrable 
[26l[2a[28l|47l|48l|49l[72l[73]. 



4 Some qualitative remarks 

It is well known that for typical complex objects like molecules, the structure 
of Raman spectra depends strongly on the mutual positions and splittings of 
excited energy levels of internal motion. In molecules a typical picture is as 
follows, cf. some pictures below based on [i^lSSlI^ : 

(i) The main background is created by the system of electronic energy levels. 
Usually they are analysed and approximately calculated on the basis of the 
Born-Oppenheimer approximation [42l [66l [67] . In principle the separation 
of those levels is such that the corresponding quantum transitions result 
in radiation of visible light. 

(m) Those levels matter of fact, bands consisting of systems of vibra- 

tional (deformative) energy levels. The frequencies of radiative transitions 
within those bands are placed within the visible light and near infrared 
ranges. 
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(Hi) And finally, the vibrational levels split into the rotational ones. Here the 
resulting frequencies belong to the far infrared and radio ranges. 





Fig. 4 

This situation is often faced with but there are exceptions, when separations 
between energy levels of various types are comparable. Then one has to do with 
some resonance phenomena known as the Jahn- Teller effect [1^ [57] . 
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Jahn-Teller effect 




Fig. 5 

And even the very Born-Oppenheimer method may be non-apphcable in such 
exceptional situations. 

In our model, when applied to molecules, fuUerens and similar objects, the 
general picture and structure of energy levels and their splittings is a bit more 
complicated and not yet analysed sufficiently in quantitative and qualitative 
details. The splitting of electronic levels into vibrational and rotational ones 
has some additional peculiarities. Namely, spatial rotations are controlled by 
the quantum number s, but quantum deformations are described by two things: 
the spectrum of and the quantum number j controlling the rotation of 
squeezing plane of the deformed object. 

Moreover, it would be difficult to estimate the structure of splittings in nu- 
clear dynamics, where, nevertheless, some interesting and nontrivial applications 
are expected. 

In our model based on the high affine symmetry one may hope that some 
partial results may be explicitly obtained. Situation certainly will be much more 
difficult for more general models of lower dynamical symmetry. The next, more 
difficult step will follow when even on the level of kinematics we give up the 
affine model of degrees of freedom and more complicated deformation modes 
are admitted. 

5 Non-afRne modes — some general and rough 
comments 

When dealing with molecular dynamics, it is quite natural to expect that affine 
modes of motion, i.e., rotations and homogeneous deformations, are most rele- 
vant for internal phenomena. They are also important in nuclear dynamics. Of 
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course, on both levels the quantization procedure must be carried out. 

Nevertheless, any object consisting of more than four material points ((n+1) 
in n-dimensional space) has also other degrees of freedom, even if in a given 
class of phenomena they are not very important and play a secondary role. In 
molecular or nuclear dynamics (and also in some macroscopic phenomena) it 
is often reasonable to establish some hierarchy of degrees of freedom starting 
from affine ones and then admitting ones more and more complicated. Degrees 
of freedom of affine motion are represented by the formula ((1]) which expresses 
the Cartesian Euler (current) coordinates as first-order polynomials of Carte- 
sian Lagrange (material) coordinates. It is natural to describe other degrees of 
freedom in such a way that Eulcr coordinates are higher-order polynomials of 
Lagrange variables [74] : 

y'{t,a) = oq'{t) + iq'A{t)a'^ + 2q\B{t)a^a^ + ■■■ 

k 

= ^pg'A,...A,(i)a^^---a-^''. (248) 

p=0 

The coefficients pQ^ AB---z(t) are generalized coordinates; one must remember 
however that they are symmetric in material (capital) indices, so to be more 
precise, independent generalized coordinates correspond, e.g., to A < B < ■ ■ ■ < 
Z . To avoid this redundancy, or rather to reduce it, it is convenient to use the 
representation in terms of the radial variable and spherical functions [74j : 

y' (^'«) = E E q'lrrdtMY'"- (a) , (249) 
1=0 m=-l \W/ 

where, obviously, |a| is the length of the material radius- vector and 

¥i^ = q'i-m, (250) 

because ?/* are real quantities. Then independent generalized coordinates cor- 
respond to TO = 0, 1, . . . , Z. If fc > 1, then on the surface of the body more than 
one deformative waves are formed. 
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k=l k>l 



Fig. 6 

This procedure is used in so-called method of virial coefRcients (widely used 
in astrophysics, cf. for instance the book by Chandrasekhar) . It was also applied 
in nuclear dynamics, where of course the quantized version of theory must be 
used [H]. 

For continuous bodies, k may be in principle arbitrary (even infinite, then 
are simply expressed as analytic functions of Lagrange coordinates). Obviously, 
for finite systems of material points, k must be finite, because one deals then 
with a finite number of degrees of freedom. The higher p oi I , the less collective 
character of the corresponding degrees of freedom, although, of course, all of 
them are "collective" in comparison to individual one-particle positions. 

Substituting ()248p into the formula of the kinetic energy, one expresses it 
through generalized velocities pQ^ AB---zit) or g^™; later on Legendre transfor- 
mation is performed to reformulate everything in phase-space terms, and finally 
the model is subject to the Schrodinger quantization. Obviously, when using 
polynomial of the order fc > 1 we loose the nice group-theoretical interpretation; 
it survives only on the level of affine (fc = 1) background phenomena. 

Remark: if fc > 1, then qq^ are not coordinates of the centre of mass any 
longer. Overlooking this fact may lead to serious mistakes. 
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